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Abstract 


We  describe  a  computer  simulation  of  atmospheric  and  target  effects  on 
the  accuracy  of  range  measurements  using  pulsed  laser  radars  (ladar)  with 
PIN  or  avalanche  photodiodes  for  direct  detection.  The  computer  simula¬ 
tion  produces  simulated  range  images  as  a  function  of  a  wide  variety  of 
environmental,  target,  and  sensor  parameters  for  ladar  with  range  accura¬ 
cies  smaller  than  the  pulse  width.  The  simulation  allows  arbitrary  target 
geometries,  and  simulates  speckle,  turbulence,  and  near-  and  far-field  ef¬ 
fects.  We  compare  simulation  results  to  actual  range  error  data  collected  in 
field  tests. 
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1.  Introduction 


Laser  radar  (ladar)  is  being  used  increasingly  to  locate  and  recognize  ob¬ 
jects  in  both  industrial  and  military  applications.  Recognition  algorithms 
often  rely  on  the  range  measurements  produced  by  ladar  rather  than 
on  intensity  measurements,  because  3-D  range  imagery  usually  contains 
more  information  than  the  intensity  images,  which  typically  contain  large 
amounts  of  speckle  and  noise.  Being  able  to  estimate  the  range  error  of  a 
ladar  configuration  is  desirable  to  help  choose  the  design  parameters  of  a 
ladar  and  to  help  design  the  recognition  algorithm.  A  number  of  papers 
address  the  problem  of  estimating  the  probability  of  detection  of  laser  re¬ 
turns  and  the  probability  density  functions  (pdf's)  of  range  error  based  on 
signal-to-noise  ratio  or  carrier-to-noise  ratio  [1],  but  these  papers  generally 
assume  a  simple  target  geometry  (e.g.,  a  plane  perpendicular  to  the  line  of 
sight)  that  allows  analytic  expressions  to  be  derived.  The  formulas  for  arbi¬ 
trarily  complex  target  geometries  cannot  be  solved  analytically,  but  require 
a  numerical  solution  that  would  be  more  computationally  complex  than 
the  simulation  described  in  this  paper.  Quantitative  analysis  is  complicated 
by  nonplanar  target  geometry,  laser  speckle,  and  atmospheric  scintillation 
effects,  requiring  a  simulation  to  be  used  instead  of  analytically  derived 
expressions. 

The  computer  simulation  that  we  describe  herein  produces  simulated  range 
images  for  a  direct  detection  pulse  ladar  under  a  wide  variety  of  conditions. 
The  simulation  incorporates  established  theory  at  each  stage.  The  simu¬ 
lation  process  is  divided  into  six  stages:  component  disassembly,  speckle 
modulation,  scintillation  modulation,  signal  assembly,  receiver  noise,  and 
pulse  detection. 

We  conducted  a  field  test  to  collect  data  to  compare  with  the  simulation  re¬ 
sults.  In  the  field  test,  we  used  two  flat  plywood  panels  painted  care  green 
as  targets.  The  range  measurements  we  collected  using  these  targets  were 
compared  to  simulation  results  to  assess  the  accuracy  of  the  simulation.  We 
intentionally  chose  a  simple  target  geometry  for  the  field  test  so  that  ana¬ 
lytical  methods  could  also  be  used  to  predict  performance  for  comparison 
to  the  simulation  results. 
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2.  Laser  Simulation 


2.1  Component  Disassembly 

To  simulate  the  effects  of  target  shape,  speckle,  and  scintillation,  we  decom¬ 
posed  the  laser  pulse  into  components  along  the  range  and  the  two  cross¬ 
range  dimensions.  Thus  the  pulse  can  be  expressed  as  U ( x ,  y,  z ),  where  x 
and  y  are  the  horizontal  and  vertical  crossrange  dimensions,  and  z  is  the 
range  dimension  (the  direction  that  the  pulse  is  traveling).  The  beam's  ir- 
radiance  profile  is  assumed  to  be  Gaussian,  and  the  pulse  shape  in  the  z 
dimension  is  an  input  to  the  model  and  can  be  written  as  V(z),  or  alterna¬ 
tively  as  V ( ct ),  where  c  is  the  speed  of  light  and  t  is  time.  So  we  decom¬ 
posed  the  outgoing  pulse  as 

U(x,  y,  z)  =  PsG(x,  y)V(z),  (1) 

where  Ps  is  the  total  pulse  power;  x,  y,  and  z  take  on  discrete  values;  G(x,  y) 
is  the  proportion  of  energy  within  a  component  located  at  (x,  y)  under  the 
two-dimensional  Gaussian  curve  at  (x,y);  and  V(z)  is  the  discrete  pulse 
shape  in  the  range  dimension  shown  in  figure  1.  The  integrals  of  G(x,  y) 
and  V (z)  are  both  one  because  they  are  normalized  to  unity. 

Each  crossrange  component  (x,  y)  corresponds  to  the  energy  in  a  5-  x  5- 
cm  square  area  (this  size  can  be  set  arbitrarily,  of  course)  at  target  range; 
thus,  the  number  of  crossrange  components  that  the  simulation  uses  de¬ 
pends  on  the  range  and  the  beam  divergence  of  the  sensor.  In  table  1  we 
show  the  portion  of  energy  in  each  crossrange  component  for  a  pulse  that 
spreads  to  an  area  of  25  x  25  cm  at  target  range.  The  distance  that  each 
component  travels  is  determined  by  a  geometric  model,  which  is  an  input 
to  the  simulation,  and  has  resolution  of  5  x  5  cm,  matching  the  crossrange 
decomposition  of  the  pulse.  At  the  target,  each  component  is  treated  as  if 
it  encounters  a  resolved  planar  surface  perpendicular  to  the  line  of  sight  to 
the  sensor.  However,  a  pulse  may  see  nonperpendicular  surfaces  and  un¬ 
resolved  surfaces  as  a  set  of  components  at  differing  ranges.  The  pulses  are 
not  split  by  frequency  because  we  assume  a  quasi-monochromatic  source. 
For  a  broadband  source,  such  as  a  semiconductor  laser,  accurate  model¬ 
ing  of  speckle  would  require  either  breaking  up  the  pulses  by  frequency, 
or  using  a  different  distribution  than  the  one  we  use  to  model  speckle  and 
turbulence  modulation. 
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Figure  1.  Simulated  laser  pulse  shape. 


Table  1.  Portion  of  total  pulse  energy  in  crossrange  components  with  Gaussian  irradiance 
profile  (G(x,  y)). 


Location 

X  =  -0.1 

x  =  -0.05 

x  =  0.0 

x  =  0.05 

*  =  0.1 

y  =  0.1 

0.003 

0.014 

0.022 

0.014 

0.003 

y  =  0.05 

0.014 

0.059 

0.095 

0.059 

0.014 

O 

O 

II 

0.022 

0.095 

0.154 

0.095 

0.022 

y  =  -0.05 

0.014 

0.059 

0.095 

0.059 

0.014 

t-H 

O 

1 

II 

0.003 

0.014 

0.022 

0.014 

0.003 

The  standard  radar  equation  governs  the  return  power  received  for  each 
component  sent  [1],  Therefore,  the  received  component  power  is  given  by 

PR  =  Pce-2aR^ep,  (2) 

where  Pr  is  the  received  component  power,  Pc  is  the  transmitted  compo¬ 
nent  power,  a  is  the  atmospheric  extinction  coefficient,  R  is  the  range  to 
the  target,  d  is  the  effective  diameter  of  the  receiver's  clear  aperture,  e  is  the 
receiver's  optical  efficiency,  and  p  is  the  diffuse  reflectivity  of  the  resolved 
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target.  The  resolved  target  assumption  is  applied  only  at  the  level  of  the 
component;  the  pulse  itself  may  be  unresolved. 

Speckle  and  Turbulence  Modulation 

The  power  received  from  each  component  is  modulated  by  speckle  and 
turbulence.  Speckle  is  applied  to  each  component  with  the  use  of  the  expo¬ 
nential  distribution  to  modulate  the  received  power,  because  we  assume 
a  quasi-monochromatic  source.  For  broadband  sources,  the  method  de¬ 
scribed  by  Parry  [2]  could  be  used  to  compute  the  pdf.  Thus,  the  power 
received  from  a  component  in  the  presence  of  speckle  is 


P speckle  P{,Pr)i 


(3) 


where  S  is  an  exponentially  distributed  random  variable  with  parameter 

Pr, 

—  s 

Prob[S(X)  =  s]  =  ^~,  (4) 

where  A  is  the  mean  of  the  exponential  random  variable  [3].  Similarly,  tur¬ 
bulence  is  applied  by  modulating  the  power  of  each  component  with  a 
lognormal  random  variable.  Some  authors  have  argued  that  other  distri¬ 
butions  for  turbulence  are  more  appropriate,  especially  in  high  turbulence 
in  which  the  K  distribution  matches  the  data  well  [4],  although  it  tends  to 
underestimate  probabilities  of  high  irradiances  [5].  Experimental  evidence 
suggests  that  in  the  presence  of  significant  aperture  averaging,  the  statis¬ 
tics  of  the  irradiance  are  lognormal  even  in  the  high-fluctuation  regime  [6]. 
Many  other  pdf's  have  been  suggested,  but  none  have  been  shown  to  fit 
the  data  under  all  conditions  [5].  The  simulation  uses  modular  code,  mak¬ 
ing  the  pdf  easy  to  change  if  consensus  is  reached.  The  mean  /t  of  the  vari¬ 
able  is  the  power  before  the  application  of  the  turbulence.  We  determine 
the  mean  normalized  variance  erf  by  [7] 


aI  —  1  target^  I  point’ 


(5) 


where  the  term  'ytarget  accounts  for  the  target  averaging  of  turbulence,  cr]point 
is  the  intensity  fluctuation  for  a  point  target  without  aperture  averaging, 
and 

alpoint  =  1.23cl(k)7H2R)%,  (6) 

where  R  is  the  range  from  the  target  to  the  sensor,  k  is  the  wave  number, 
and  (?n  is  the  turbulence  refractive  index  structure  constant,  which  is  depen¬ 
dent  on  atmospheric  conditions.  The  value  of  (?n  used  in  the  simulation  is  a 


user  controlled  input  parameter.  We  obtained  typical  values  from  Shapiro 
et  al  [1]  and  determined  the  target  averaging  term  7 target  by  [7] 


O' target 


(7) 


where  the  value 


Teff  —  min  \  ftgt  1  f'beam  (P)  :■  Pfov  (  7?)  ] 


(8) 


is  a  measure  of  the  effective  averaging  area  and  pi  is  the  long-term  turbu¬ 
lence  cell  size.  We  calculated  the  pi  using 


Pi  = 


(9) 


where  [7] 


P  o  -  P  o 


(s) 

Po 

( v ) 

Po 


ZB 

fxmt 

5 


(p) 


P-&l2  +  (&2(1  +  lT>r&  _ 
Ll-¥(^)  +  ¥(7t)2  +  i(£)2^J 

3 


3  Vmt 


(0.5A:2c^i?)“6“, 
(1.46/c2c2i?)“5“, 
krbo 


=  (^) 


1 


rf1)2 

—^60 


+ 


4rbJ 


tan(0bl/2)-tan(^_)’ 
2 TfeO 

(p)  ' 

Po 


(10) 

(11) 

(12) 

(13) 

(14) 

(15) 


The  Rytov  solution  (eq  5)  predicts  that  the  variance  of  the  intensity  fluctu¬ 
ation  increases  indefinitely  as  range  or  the  structure  constant  c2  increases. 
Empirical  measurements  show  that  the  normalized  variance  of  intensity 
fluctuations  saturate  at  approximately  1,  whereas  the  normalized  irradi- 
ance  variance  is  unbounded  in  the  Rytov  solution  [8].  The  empirical  curve 
of  this  saturation  is  shown  in  figure  2.  This  curve  is  stored  as  a  lookup  table 
in  the  simulation. 

The  received  power  in  a  component  after  turbulence  is  then 


Pt 


turbulence 


—  eG{a,b)  p 


speckle  i 


(16) 
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Figure  2.  Saturation  of  intensity  fluctuation. 

where  G  is  a  Gaussian  random  variable  with  mean  a  and  variance  b.  We 
determined  the  parameters  a  and  b  from  the  mean  and  variance  of  the  log¬ 
normal  distribution  using 


b 

a 


I  log( 


a- 


pe2  +  l'’ 


log(M)  -  y- 


(17) 

(18) 


The  speckle  and  turbulence  modulation  of  the  component  intensities  is  not 
independent  from  component  to  component.  The  correlation  length  of  the 
speckle  and  turbulence  in  the  crossrange  dimension  are  calculated,  and 
components  are  grouped  together  so  that  the  combined  components  are 
the  size  of  a  turbulence  or  speckle  cell.  The  turbulence  cell  size  is  equiva¬ 
lent  to  the  value  pi  calculated  in  equation  (9).  The  speckle  cell  size  in  the 
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target  plane  is 


l speckle  —  r-r—  i  (19) 

V  A-o 

where  Aa  is  the  area  of  the  receiver  optics.  The  correlation  time  of  speckle 
and  turbulence  for  a  single  path  is  assumed  to  be  greater  than  the  pulse 
width,  so  that  turbulence  and  speckle  for  a  given  pulse  is  a  function  only 
of  crossrange  dimensions  (x,  y),  not  of  the  range  dimension  z. 

Once  we  applied  speckle  and  turbulence,  we  summed  the  components  in 
the  crossrange  dimension  to  form  the  return  signal  at  the  detector.  Each 
component  is  shifted  in  range  corresponding  to  the  range  to  the  target  for 
that  component.  The  returning  power  at  the  detector  then  has  the  form 

Ur(z)  =  J2Pturbulence(x,y)V(z  -  2 T(x,y)),  (20) 

x,y 

where  the  crossrange  dependence  of  Pt,UrbuiencAx,  y )  is  made  explicit,  and 
T(x,y)  represents  the  range  to  the  target  at  location  (x,y). 


2.3  Detector  and  Amplifier  Noise 

The  return  signal  power  is  multiplied  by  (arm,  which  is  a  function  of  the 
geometric  radius  of  the  laser  spot  image,  turbulence  blur  circle,  and  diffrac¬ 
tion  limited  blur  radius.  We  then  calculated  the  factor  as  [9] 


€a(R,r) 
A  (a,  b,  c) 
^(b,  a,  c) 
Fgc(R) 


r'Wjr*. 

vr r{;{R)  Jo 

A  ( r  raper  Jrcvr  V  frcvr  \  A  L  Fpbsfrcvr  r  f  \ 

IYvDi  R  ’  R  )  R  >  R) 

aPer  Jrcvr  ^  2 


R 

62'E(a,  6,  c)  +  a2'V(b.  a,  c )  —  ac sin('k(6,  a,  c)), 
'C2  +  b2  —  a2 


arccos 


nun 


/c  +  o  —  a  \ 

V  2 be  A 

(  — — l) 

^Ablur(R)’  )' 


(21) 

(22) 

(23) 

(24) 

(25) 


where  Ad  is  the  area  of  the  detector,  raper  is  the  radius  of  the  receiver  aper¬ 
ture,  r0bs  is  the  radius  of  the  obscuration,  Fgc(R)  is  the  detector  geometric 
compression  factor,  and 


Ablur(Ffy  (-^)  •> 

(26) 

^blur^R)  ^peom(-R)  T'turb+dif /(-^) * 

(27) 
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The  diffraction  limited  blur  radius  and  turbulence  blur  circle  radius  for  the 
short  exposure  case  can  be  taken  as 


Vturb+dif  f 


rdiff 


_  „  [l  I  rdiff  2]  2 

-  rdiff  1  +  j 

(28) 

L  rtrb  J 

_  1.22  A  frcvrQ 

drcv 

(29) 

2 

T'turb  i 

(30) 

TTlStrb 


where  vtrb  is  the  turbulence  cutoff  spatial  frequency  given  by  the  solution 
v  of  the  equation  [10] 


( 


1 


)!(- 


3.44  A/, 


)  =  v  1  —  a( 


A/, 


(31) 


and 

r0  =  2.1(1.46£:2c;;i?)-f  (32) 

is  Fried's  coherent  aperture  diameter  because  of  turbulence.  The  geometric 
radius  of  the  laser  spot  image  is 


fr, 


-n(R), 


'  geom.  ~  ^ 

where  the  laser  spot  radius  rb(R)  is  calculated  as  [11] 


(33) 


rb(R)  = 


R2 

k2rbo 


2  ,  R  ,  2  4 R 2 

+  r“<1_S^)  +W,\ 


if  R  >  A;min(/^,  4rf0)  or  pi  >  2rbo,  and 


rb(R )  = 


R2 


*2  r 


k2- 2 


bO 


2  .  R  ,2  4  R 

+  r6o(i-  —  7  )  +  T2“2 

Jxmt  P[ 


1  -  o 


m(—) 

V2rw/ 


(34) 


(35) 


fxrnt  =  - — - y - ,  (36) 

tan(V2)-tan(^)’ 

otherwise.  In  equations  (31)  to  (36),  rb o  is  the  radius  of  the  beam  waist, 
frcvr  is  the  focal  length  of  the  receiver,  fxrnt  is  the  effective  focal  length  of 
the  transmitter,  0hl/2  is  the  half-angle  transmitted  beam  divergence,  Q  is 
the  quality  factor  of  the  optics,  drcv  is  the  diameter  of  the  effective  clear 
aperture,  A  is  the  wavelength  of  the  laser,  k  is  the  wave  number,  and  R  is 
the  range  to  the  target. 
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The  return  signal  is  subject  to  background  noise,  shot  noise,  amplifier  noise, 
and  dark  current  noise,  all  of  which  are  independent,  identically  distributed 
(iid)  Gaussian.  The  detector  noise,  consisting  of  dark,  shot,  and  background 
noise,  is  calculated  as  [12] 


A  /'.  1  ’itil  iflay 


\/2  ecB[Ids  +  (Idb  +  Ib  +  IS)MW\ 

n 


(37) 


where  1Z  is  the  responsivity  of  the  detector,  B  is  the  electrical  bandwidth. 
Ids  is  the  surface  dark  current.  Bib  is  the  bulk  dark  current  at  unity  gain, 
lb  =  IZPb  is  the  current  because  of  the  background  illumination,  Is  = 
IZPturbuience  is  the  current  because  of  the  received  signal,  M  is  the  detec¬ 
tor  gain,  and  F  is  the  excess  noise  factor  because  of  the  detector  gain.  Pb  is 
the  background  power  calculated  from 


Pb 


pi  s  u  n  If  ■'  1  r  S 1  n 


2 

AA, 


(38) 


where  Ar  is  the  area  of  the  receiver,  Tr  is  the  transmission  of  the  receiver,  p  is 
the  background  reflectance,  6 fov  is  the  field  of  view,  hsun  is  the  background 
solar  irradiance,  and  AA  is  the  optical  bandwidth. 

Amplifier  noise  is  calculated  assuming  a  basic  RC  filtered  amp.  It  follows 


NEPa„ip  = 


4  kTBN 
RlP2 


(39) 


where  k  is  Boltzmann's  constant,  T  is  the  temperature  in  Kelvins,  B  is  the 
electrical  bandwidth,  N  is  the  noise  factor  for  the  electronics,  and  Rl  is  the 
load  resistor  calculated  from 


Rl 


1 

2  irBC’ 


(40) 


where  C  is  the  capacitance  of  the  detector. 

The  noises  are  added  in  quadrature  to  determine  the  overall  noise  figure 
for  the  detectors 


NEPTotaI  =  jNEP\mr  +  NEPiFB  •  («> 

The  NEP  is  the  standard  deviation  of  the  Gaussian  distribution  of  the  ad¬ 
ditive  noise  in  the  optical  power  domain.  The  simulation  generates  a  Gaus¬ 
sian  random  variable  with  zero  mean  and  a  standard  deviation  equal  to 
NEP.  This  Gaussian  noise  is  added  to  the  signal,  and  then  is  filtered  ac¬ 
cording  to  the  electrical  bandwidth  of  the  sensor.  The  noisy  signal  is  then 
passed  to  the  pulse  detector. 
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The  pulse  detection  portion  of  the  simulation  is  modifed  according  to  the 
sensor  being  simulated.  Commonly  used  pulse  detection  algorithms  in¬ 
clude  matched  filtering  followed  by  peak  detection,  threshold  detection  of 
the  rising  edge  of  the  pulse,  or  detection  averaging  of  the  rising  and  falling 
edge  of  the  pulse.  Quantization  error  of  the  system  clock  is  simulated  by 
shifting  the  transmitted  and  received  pulse  by  a  random  variable  that  is 
uniformly  distributed  across  ±  one-half  of  a  clock  cycle. 

Figure  3  shows  the  effect  of  the  noise  on  the  simulated  received  pulse.  Fig¬ 
ure  4  shows  a  simulated  ladar  range  image  exhibiting  the  effects  of  the 
signal  fluctuations  and  noise  included  in  the  simulation. 
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Figure  4.  Input  range  image  and  resulting  ladar  image.  (Black  pixels  indicate  dropouts  and 
white  pixels  indicate  anomalies.) 
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3.  Comparison  of  Simulated  Data  to  Real  Data 


As  a  check  on  the  quality  of  the  simulation,  we  gathered  a  small  amount  of 
real  data.  We  painted  a  4-  x  8-ft  sheet  of  plywood  care  green  and  positioned 
it  perpendicular  to  the  line  of  sight  to  the  sensor.  We  captured  ladar  images 
of  the  plywood  and  extracted  an  empirical  pdf  from  the  images.  Figure  5 
shows  simulated  range  error  as  a  function  of  range.  The  two  x's  mark  the 
range's  standard  deviation  for  the  real  ladar.  Additional  simulation  runs 
showed  that  the  ladar  is  clock  quantization  error  limited  at  short  ranges  and 
signal-to-noise  ratio  limited  at  longer  ranges.  Many  ladar  systems,  includ¬ 
ing  the  Lockheed  Martin  Vought  system,  avoid  the  sharp  increase  in  error 
at  longer  ranges  by  returning  no  range  value  for  return  pulses  that  do  not 
exceed  several  times  the  ambient  root  mean  square  (rms)  noise  level.  Such 
a  return  pulse  is  called  a  dropout.  The  percentage  of  dropout  pixels  would 
then  increase  as  range  increases,  while  the  range  error  on  nondropout  pix¬ 
els  would  increase  much  more  slowly.  Table  2  gives  a  partial  list  of  sensor 
parameters. 

The  plywood  target  data  that  we  collected  have  characteristics  that  limit 
validation  uses.  The  Lockheed  Martin  Vought  ladar  has  a  relatively  large 
aperture,  so  aperture  averaging  makes  the  effects  of  speckle  and  scintil¬ 
lation  almost  negligible.  The  detector  noises  are  dominated  by  amplifier 
noise,  so  background  and  shot  noise  effects  are  not  adequately  tested.  While 
the  comparison  of  real  and  simulated  data  was  useful  in  validating  the  sim¬ 
ulation,  it  is  far  from  thorough.  A  thorough  empirical  validation  would  re¬ 
quire  laser  systems  having  at  least  several  different  values  of  each  sensor 
parameter,  so  that  each  sensor  parameter  is  accurately  modeled. 

To  test  the  speckle  and  scintillation  portion  of  the  simulation,  we  had  to 
reduce  the  aperture  size  until  aperture  averaging  was  insufficient  to  make 
speckle  and  scintillation  effects  negligible.  The  aperture  size  was  reduced 
while  the  output  power  of  the  laser  was  increased  so  that  average  return 
power  remained  a  constant;  thus  the  primary  cause  of  variation  in  perform¬ 
ance  is  the  reduction  in  aperture  averaging  that  occurs  at  smaller  aperture 
sizes.  Figure  6  shows  the  resulting  simulated  performance. 
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Figure  5.  Simulated  range  error  versus  range  (x's  mark  real  data  points). 


Table  2.  Sensor  description. 


Parameter 

Value 

Laser  type 

Nd:YAG 

Wavelength 

1.574  /im 

Beam  divergence 

125  /rrad 

Pulse  width 

10  ns 

Pulse  energy 

250  /rj 
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4.  Conclusions 


We  have  described  a  computer  simulation  that  estimates  ladar  range  per¬ 
formance  as  a  function  of  sensor,  environmental,  and  target  parameters.  By 
breaking  the  pulse  into  components,  the  simulation  is  able  to  simulate  laser 
returns  for  arbitrary  target  geometries. 

Comparison  of  simulated  data  to  a  limited  set  of  real  data  shows  a  close 
match.  Unfortunately,  we  were  unable  to  collect  a  wider  variety  of  data, 
such  as  data  at  longer  ranges.  Additional  data  validation  is  needed. 

The  simulation  has  been  used  in  performance  versus  cost  trade-off  stud¬ 
ies  during  the  design  phase  of  the  Communications-Electronics  Command 
(CECOM)  Multifunction  Laser  System  (MFLS).  Currently,  it  is  being  used 
to  estimate  range  error  as  an  aid  to  designing  ladar  automatic  target  rec¬ 
ognizer  algorithms  in  a  joint  program  at  CECOM  and  the  Army  Research 
Laboratory  (ARL). 
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